Rayleigh Imaging of Graphene and Graphene Layers 
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We investigate graphene and graphene layers on different substrates by monochromatic and white- 
light confocal Rayleigh scattering microscopy. The image contrast depends sensitively on the dielec- 
tric properties of the sample as well as the substrate geometry and can be described quantitatively 
using the complex refractive index of bulk graphite. For few layers (<6) the monochromatic contrast 
increases linearly with thickness: the samples behave as a superposition of single sheets which act as 
independent two dimensional electron gases. Thus, Rayleigh imaging is a general, simple and quick 
tool to identify graphene layers, that is readily combined with Raman scattering, which provides 
structural identification. 
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Graphene is the prototype two dimensional carbon 
system [J. Its electron transport is described by the 
(relativistic-like) Dirac equation and this allows access 
to the rich and subtle physics of quantum electrodynam- 
ics in a relatively simple condensed matter experiment 
0, i, i, i, 0. The scalability of graph ene devices to 
true nanometer dimensions [?, 8, "9] makes it a promising 
candidate for future electronics, because of its ballistic 
transport at room temperature combined with chemical 
and mechanical stability. Remarkable properties extend 
to bi-layer and few-layers [5, 6, l(?,fTl|,lT2l. More funda- 
mentally, the various forms of graphite, nanotubes, buck- 
y balls can all be viewed as derivatives of graphene. 

Graphene samples can be obtained from micro- 
mechanical cleavage of graphite [2]. Alternative proce- 
dures include chemical exfoliation of graphite 
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J or epitaxial growth by thermal decomposition 

of SiC 0, H, [l9, 20]. The latter has the potential of 
producing large-area lithography compatible films, but is 
substrate limited. It is hoped that in the near future effi- 
cient large area, substrate independent, growth methods 
will be developed, as it is now the case for nanotubes. 

Despite the wide use of the micro- mechanical cleavage, 
the identification and counting of graphene layers is still a 
major hurdle. Monolayers are a great minority amongst 
accompanying thicker flakes [1]. They cannot be seen 
in an optical microscope on most substrates. Currently, 
optically visible graphene layers are obtained by placing 
them on the top of oxidized Si substrates with typically 
300 nm Si02. Atomic Force Microscopy (AFM) is viable 
but has a very low throughput. Moreover, the different 
interaction forces between the AFM probe, graphene and 
the Si02 substrate, lead to an apparent thickness of 0.5- 
1 nm even for a single layer [2, 5], much bigger of what 
expected from the interlayer graphite spacing. Thus, in 
practice, it is only possible to distinguish between one 
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FIG. 1: Schematic experimental set-up for combined Rayleigh 
and Raman spectroscopy. The inset shows a cross sectional 
view of the interaction between the optical field and graphene 
deposited on Si covered with Si02. 



and two layers by AFM if graphene films contain folds 
or wrinkles [l, [HI . High resolution transmission electron 
microscopy is the most direct identification tool [^T, ^25] , 
however, it is destructive and very time consuming, being 
viable only for fundamental studies (2l[. 

Optical detection relying on light scattering is espe- 
cially attractive because it can be fast, sensitive and not- 
destructive. Light interaction with matter can be elastic 
or inelastic, and this corresponds to Rayleigh and Raman 
scattering, respectively. Raman scattering has recently 
emerged as a viable, non-destructive technique for the 
identification of graphene and its doping [22, 23]. How- 
ever, Raman scattered photons are a minority compared 
to those elastically scattered. Here we show that the 
elastically scattered photons provide another very effi- 
cient and quick means to identify single and multi-layer 
samples and a direct probe of their dielectric constant. 

Rayleigh scattering was previously used to monitor 
size, shape, concentration and optical properties of nano- 
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particles, carbon nanotubes and viruses [2J, |25|, |26|, |27|, 
[28]. Rayleigh scattering experiments can be performed 
using two different strategies. In one, the background 
signal is minimized by making free-standing samples, as 
done in the case of carbon nanotubes [13, HJ, or by dark- 
field configurations [l^. Alternatively, the background 
intensity is utilized as a reference beam, while the saniple 
signal is detected interferometrically |23, [H, 0, [sO, [lH . 
Here, we combine the second approach with the inter- 
ferometric modulation of the contributing fields and we 
show that the presence of a background is essential to 
enhance the detection of graphene over a certain wave- 
length range. 

Graphene samples are produced by micro-mechanical 
cleavage of bulk graphite and deposited on a Si substrate 
covered with 300 nm Si02 (IDB Technologies LTD). The 
sample thickness is independently confirmed by a combi- 
nation of AFM and Raman spectroscopy. AFM is per- 
formed in tapping mode under ambient conditions. Ra- 
man spectra are measured at 514 nm using a Renishaw 
micro- Raman 1000 spectrometer. Rayleigh scattering is 
performed with an inverted confocal microscope, FiglH 
Either a He-Ne laser (633 nm) or a collimated white-light 
beam are used as excitation source. Coherent white-light 
pulses are generated by pumping a photonic crystal fi- 
bre with the output of a Ti:Sa oscillator operating at 760 
nm. The beam is refiected by a beam splitter and focused 
by a microscope objective with high numerical aperture 
(NA= 0.95). However, the objective lens is not totally 
filled, which results in an effective NA~0.7 thereby in- 
creasing the image contrast as discussed at the end of this 
paper. The scattered light from the sample is collected in 
backscattering geometry, transmitted by a beam splitter 
and detected by a photon-counting avalanche photodiode 
(APD), FigUl Alternatively, the refiected light is filtered 
using a notch filter to remove the laser excitation and sent 
to a spectrometer. This allows simultaneous Rayleigh 
and Raman measurements^ Fig jl|2k . Confocal Rayleigh 
images are obtained by raster scanning the sample with 
a piezoelectric scan stage. The acquisition time per pixel 
varies from few ms in the case of Rayleigh scattering to 
few minutes for Raman scattering. This empirically in- 
dicates that Rayleigh measurements are almost 5 orders 
of magnitude quicker than Raman measurements. The 
spatial resolution is ~ 800 nm. 

FigEl^b) shows an AFM image of monolayer graphene. 
The AFM cross section gives an apparent height of ~0.6 
nm. Raman spectroscopy confirms that the sample is a 
single layer (FiglSfa)) [221. FiglSfb) is the correspond- 
ing confocal Rayleigh image obtained with monochro- 
matic laser light (633 nm). Fig|3l[a) shows an optical 
micrograph of a sample composed of a varying number 
of layers. Once the single layer is identified by Raman 
scattering, we get the total number of layers from the 
measured AFM height, considering the interlayer spac- 
ing of ~ 0.33 nm: z [nm]= 0.27 + 0.33 N. This confirms 




1600 2000 2400 2800 

(b) Raman shift (cm"^) 



5 |im 



(c) 



# 

S.OLim 



FIG. 2: (a) Raman spectrum at 514 nm, showing the features 
of graphene (b) AFM image of single layer graphene 

(c) Confocal Rayleigh image obtained by raster scanning the 
sample with a piezoelectric scan stage. 



that the sample is composed of 1, 2, 3 and 6 layers, as 
for Figl3] (a). These layers have a slightly different color 
in the optical microscope (Figl3](a)). It appears that the 
darker color corresponds to the thicker sample. Note, 
however, that the color of much thicker layers (more than 
10 layers) does not follow this trend and can change from 
blue, to yellow, to grey. The number of layers is further 
confirmed by the evolution of the 514 nm Raman spec- 
tra [22], Fig. [3] (b). FiglH^a) shows a confocal Rayleigh 
map for 633nm excitation. The signal intensity of in 
FigHl appears to increase with N. 

We now discuss the physical origin of the image con- 
trast {6). This is defined as the difference between sub- 
strate and sample intensity, normalized to the substrate 
intensity. The single layer contrast at 633 nm is ~ 0.08. 
The contrast is positive, i.e. the detected intensity from 
graphene is smaller than that of the substrate. The 
Rayleigh images in FigO (c) and FigH] (a) are reversed 
for convenience, in order to compare them with AFM. 

We explain the sign and scaling of the contrast for 
increasing N in terms of interference from multiple re- 
fiections. The inset in FigH] shows a schematic of the 
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FIG. 3: (a) Optical micrograph of multi- layer with 1, 2, 3 and 
6 layers; (b) Raman spectra as a function of number of layers. 
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FIG. 4: (a) Three-dimensional confocal Rayleigh map for 
monochromatic 633 nm excitation. The window size is 49/im 
X 49/im; (b) Experimental (dots) and theoretical (line) con- 
trast as a function of excitation wavelength. 



interaction between the light and graphene on Si+Si02. 
When the light impinges on a multi-layer, multiple reflec- 
tions take place [32|. Thus, the detected signal (I) results 
from the superposition of the reflected field from the air- 
graphene (£^g), graphene-Si02 (^siOs): and Si02-Si in- 
terfaces (Esi). The back-ground signal (Ibq) results from 
the superposition of the reflected field from the air-Si02 
interface and the Si substrate. 

Before giving a complete quantitative model, it is use- 
ful to consider a simplified picture that captures the 
basic physics and illustrates why a single atomic layer 
can be visualized optically. The field at the detector is 
dominated by two contributions: the reflection by the 
graphene layer, and the reflection from the Si after trans- 
mission through graphene and after passing through the 
Si02 layer twice. Thus, the intensity at the detector can 
be approximated as: 



\Eg\^ ^\Es^\^ ^2\EG\\Es^\ COS cl, (1) 



where <p is the total phase difference. This includes the 
phase change due to the optical path length of the oxide, 
dsi02i that due to the reflection at each boundary, 
i&Si and ^G' 



where nsi02 is the refractive index of the oxide and Aq 
is the wavelength of the light in vacuum. Assuming the 
field reflected from graphene to be very small, \Eg\^ ^ 0, 
the image contrast 5 results from interference with the 
strong field reflected by the silicon: 

5 = {Is^ - I)/Is^ ^ -2 • \EG\/\Es^\ ■ COS (/> (3) 

The sign of 5 depends on the sign of cos ^, which is given 
by Eq. [2l The reflectance, R, is the ratio between the 
reflected power to the incident power [sS] . Assuming the 
Si reflectance as one, Eq. [3] can be written as: 



-2^Rg cos( 



(4) 



where Rg is the reflectance of graphene. This is in turn 
related to the reflection coefficient tg [32[: 



TG = V^G • exp{ii^G) 



(5) 



- i^Si + 27r nsi022dsi02/>'o) 



(2) 



Eq. [4] shows that the main role of the Si02 is to act as a 
spacer: the contrast is defined by the phase variation of 
the light refiected by the Si [33] . Thus, the contrast for a 
given wavelength can be tailored by adjusting the spacer 
thickness or its refractive index. 

In order to investigate the wavelength dependence of 
the image contrast, we perform Rayleigh spectroscopy 
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with a white-light source. A grating is used to analyze 
the detected light. Figl4](b) shows that for N=l the con- 
trast is maximum at ~ 570 nm. The contrast at 633 nm is 
~ 0.08, in agreement with the monochromatic Rayleigh 
scattering experiment. The contrast is zero at 750 nm 
and it is small and negative for A> 750 nm. From Eqs.[2] 
and [4] and assuming i9si = — tt, the phase of graphene 
is i^G — — TT as expected for an ultra-thin film [32]. The 
contrast decreases in the near IR (for dsi02= 300 nm) 
since the wavelength becomes larger than twice the opti- 
cal path length provided by the Si02-spacer. Fig. [H (b) 
shows that while the contrast increases for increasing N, 
the phase remains constant. 

We now present a more accurate model, with no as- 
sumptions, which describes the light modulation by mul- 
tiple reflections based on the recurrent matrix method 
for reflection and transmission of multilayered films [34] . 
We calculate the total electric and magnetic fields in the 
various layers, applying the boundary conditions at every 
interface. The fields at two adjacent boundaries are de- 
scribed by a characteristic matrix. This depends on the 
complex refractive index and the thickness of the film 
and the angle of the incident light [34]. By computing 
the characteristic matrix of every layer and taking into 
account the numerical aperture of the objective and the 
filling factor, it is possible to find the reflection coefficient 
for an arbitrary configuration of spacer (2) and substrate 
(3) and for any number of graphene layers (G). Assuming 
two counter-propagating waves, the standard boundary 
conditions for the reflection coefficient of a normally in- 
cident wave is: 
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with = 27rnGdG / and ^2 = 27rn2C^2/Ao- For inci- 
dence at an angle 0, with s-polarization (transverse elec- 
tric field), the same formula applies with the substitu- 
tion Hi Hi COS Oi^ while for p-polarization every ratio 
changes Ui / Uj ni cos Oj / rij cos Oi . The phases change 
in both s and p polarizations to (j)G = ^TTUGdG cos Og / Ao 
and (j)2 = 27rn2(i2 cos^2/Ao. The angle Oi for every layer 
is obtained from Snell's law: Oi = arcsin(sin^o/^i)- In 
case any of the layers is absorbing (as in graphene and 
Si), we need use an effective index n- = firii.Oo) which 
depends on the incident angle from vacuum ^0 32, 35[. 
In this case the corresponding refraction angle is Oi = 
arcsin[sin ^0 / ^^(n^] • 

The matrix method requires as input the complex re- 
fractive index of the sample. The frequency dependent Si 
and Si02 indexes are taken from Ref. [s^. For graphene, 
few layers graphene and graphite, this is anisotropic, de- 
pending on the polarization of the incident light. For 
electric field perpendicular to the graphene c-axis (in- 
plane) we need n^Perp, while for electric field parallel 
to the c-axis we need riGParai- To get these, we use 
the experimental refractive index taken from the elec- 
tron energy loss spectroscopy measurements on graphite 
of Ref. [37]. For s-polarized light (electric field restricted 
in the plane) the refractive index to be used is simply 

= TiGPerp- For p-polarizatiou, both in-plane and out- 
of-plane field components exists. Thus we have an an- 
gle dependent refractive index = ^^pg^^ cos^ Oi + 
^Q^p^^^isw? Oi^ where the refracted angle Oi has to be 
calculated self-consistently with Snell's law. In order to 
account for the numerical aperture in the experiment, 
we need to integrate the response of all possible inci- 
dent angles and polarizations with a weight distribu- 
tion accounting for the Gaussian beam profile used in 
the experiment /(6>o) = e~^^^" 6'0/sin ^rr^27rsin^o, where 
Om = arcsin{NA). 

FiglH^b) shows the calculated contrast for N between 
1 and 6 (lines). This is in excellent agreement with the 
experiments: i) the contrast scales with number of lay- 
ers; ii) it is maximum at ~570 nm; iii) no phase shift is 
observed in this N range. Thus, for N between 1 and 6, 
cos0(A = 570 nm) = — 1. The contrast of graphene at 
570 nm is 0.1. From Eqs. H and[5] we get tg (A= 570 
nm)= 0.05. Thus, Rg (A= 570 nm)= 0.003. 

It quite remarkable that, without any adjustable pa- 
rameter, graphene's response can be successfully mod- 
eled using graphite's dielectric constant. This implies 
that the optical properties of graphite do not depend on 
the thickness, i.e. graphene and graphite have the same 
optical constants. The electrons within each graphene 
layer form a two dimensional gas, with little perturba- 
tion from the adjacent layers, thus making multi-layer 
graphene optically equivalent to a superposition of al- 
most non-interacting graphene layers. This is intuitive 
for s-polarization. However, quite notably this still holds 
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FIG. 5: (a) Maximum contrast at 633nm as a function of 
N; (b) Calculated contrast of graphene as a function of oxide 
thickness and excitation wavelength. Dotted lines trace the 
quarter-wavelength condition. 



when the out-of-plane direction (p-polarization) is con- 
sidered. This is because, compared to the in-plane case, 
graphite's response is much smaller, and in addition it 
gets smeared out by the NA integration. Thus, the max- 
imum contrast (A= 570 nm) of a N-layer is: S{N) = 
0.1 ■ N. Fig. 6(a) shows that this approximation fails for 
large N. When valid, the relation between topography 
and contrast is given by: z[nm]= 0.27 + 3.3S{N). 

FiglSJb) plots the contrast as a function of wavelength 
and Si02 thickness for a single layer. The maximum con- 
trast occurs at the minima of the background reflectivity. 
This is expected because this is the most sensitive point 
in terms of phase matching, and small changes become 
most visible. Thus, the optimal configuration requires 
the Si02 to be tuned as an anti-reflection (AR) coat- 
ing, i.e. with its optical length a quarter wavelength. 
The yellow dotted lines trace the quarter- wave condition 
2^5^02 <^S'i02/'^o = (^ + 1/2), and indeed they closely fol- 
low the calculated contrast maxima. A second point of 
interest are the bright spots around 275 nm. These are 
due to the absorption peak at the tt ^ tt* transition of 
graphite [37]. For this excitation, the graphene mono- 
layer not only becomes much more visible, but the con- 
trast change also directly reveals the frequency depen- 
dence of the graphene 's refractive index. Thus, as for 



nanotubes [27|, |28|, white light Rayleigh scattering is a 
direct probe of the dielectric function. 

For thicker samples {N > 10) the phase change due to 
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FIG. 6: (a) Calculated contrast of 50 layers as a function of 
oxide thickness and excitation wavelength; (b) contrast at 633 
nm for 300 nm Si02 as a function of N; (c) Experimental 
contrast at 633 nm for a thick sample. 



the optical path in graphite cannot be neglected. FiglHl^a) 
shows the calculated contrast for a 50 layer sample as a 
function of Si02 thickness, while FiglBJb) plots the con- 
trast for a fixed 300nm Si02 thickness, but for a variable 
number of layers. At 633 nm, as N increases, the response 
first saturates, then decreases and red-shifts, finally be- 
coming negative, as found experimentally (Fig. [6] (c)). It 
is also interesting to note that for small N the variation 
along the vertical (wavelength) axis is largely between 
zero and positive (i.e. reflectivity reduction only), while 
for large number of layers, the variation is from positive 
to negative (i.e. both reflectivity reduction and enhance- 
ment). This points to two different mechanisms. For 
small N, the effect of the graphene layers is just to change 
the reflectivity of the air/Si02 interface, while they offer 
no significant optical depth. For large N, on the other 
hand, the reflectivity of the air/graphene interface satu- 
rates while the effect of the increasing optical path within 
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the now thick graphite layer becomes significant. This 
change is not a monotonic function of N. While these 
two effects are different, they both contribute to a shift of 
the refiectivity resonance condition, and thus explain the 
increasing opaqueness of thicker graphene layers, when 
measured for a fixed excitation energy. 

It is also interesting to consider the contrast as a func- 
tion of NA. The calculations show that measurements 
at a reduced NA would give a stronger contrast, as one 
could intuitively expect. However, there is a nontrivial 
implication when varying NA, if one tries to maximize 
the contrast by using the anti-refiection coating rule for 
the spacer. The ideal AR coating over a substrate of 
index Usubst must have an index n spacer = V'^subst and 
quarter wave thickness dspacer = (m + 1 / 2)\q / 2^Jnsubst' 
Since nsi ~ 4 at 600 nm, it is natural to think that 
a spacer of n=2 (e.g. Si3N4) would be ideal. To ex- 
plore this. Fig. [3 plots the contrast for different NAs as 
a function of Uspacer at 600nm and for spacer thickness 
dspacer = ^00nm{nsi02 / 'Ti spacer) , which scrvcs to main- 
tain the AR condition and thus the maximum response. 

Contrary to expectations, the contrast maximizes for 
different spacer indexes depending on NA. For normal 
incidence, it is maximum at 1.93 with a huge contrast 
of 0.6 for a single layer. Fig. [3 It also has a strong 
variation thereafter, and becomes negative. As NA fur- 
ther increases, the peak moves to a smaller index (around 
1.5 for NA=0.7), becomes relatively fiat, and eventually 
goes to Uspacer = 1- Thus, for large NA, it makes lit- 
tle difference what the spacer index is, as long as the 
quarter- wave condition is satisfied. Indeed, for the ideal 
AR condition the background refiectivity goes to zero 
and thus the contrast becomes large, however this condi- 
tion strongly depends on the incidence angle and is thus 
easily destroyed at large NAs. For all possible spacer 
refractive indexes, a reduction in NA results into an in- 
creased contrast, however, the magnitude of this increase 
varies: at n=1.5 going from 0.7 to 0.0 NA changes the 
contrast by a factor 2, while at n=1.9 one can gain a fac- 



FIG. 8: (a) Optical micrograph of flakes on glass. (b) 
Rayleigh image at 633 nm excitation. The contrast is much 
higher compared to (a). 

tor of 6, Fig. [3 For maximum visibility, a SisN^ spacer 
of thickness 225nm with NA=0.0 would be ideal. How- 
ever, if high resolution is needed, as for nano-ribbons or, 
in general, to analyze edges and defects, a compromise 
between resolution and image contrast is necessary. 

A second point to note is that for all NAs the con- 
trast converges to the same value for n=l, i.e. for a sus- 
pended graphene layer over the substrate. Indeed, opti- 
cally visible suspended layers were recently reported (see 
Fig.l of Ref. [38]). Maximum visibility is achieved if the 
quarter- wave condition is satisfied, as indeed in Ref. [s^, 
where the 300 nm Si02 spacer is etched to create an air 
gap between graphene and the Si substrate. Interest- 
ingly, in this case any measurement with any NA will 
yield the same contrast. The same considerations are 
relevant for the case of a thin free-standing spacer (no 
substrate). By tuning at the low refiection point (now at 
half- wavelength) and with an NA=0.0 one could get fair 
contrasts. However, as soon as NA increases, the reso- 
nance condition is destroyed and the contrast becomes 
much smaller than for the Si02/Si system. 

The matrix method can be extended to every film con- 
figuration. To prove this, we measure graphene layers 
on glass. For N=l, the calculated contrast at 633 nm is 
expected to be~-0.01. Note the different sign compared 
with the Si/Si02 substrate. This is due to the differ- 
ent optical properties of glass and Si. FiglHl^a) shows 
an optical micrograph of a multi-layer and FiglHJb) the 
corresponding Rayleigh image at 633 nm. Raman spec- 
troscopy shows that the sample is composed of layers 
of different thickness: A (7-10 layers), B (3-6 layers), C 
(1-2 layers). Although the contrast is lower compared 
to Si/Si02, Rayleigh spectroscopy allows a better con- 
trast and resolution compared with the optical micro- 
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FIG. 9: Calculated contrast of graphene on glass at different 
wavelengths. 



scope, where 10 layers are already difficult to detect and 
single layers are practically invisible. Note that the use 
of UV light could enhance the contrast to ~ -0.04 at 300 
nm excitation (Fig. [9](b)). 

In conclusion, we used white light illumination com- 
bined with interferometric detection to study the contrast 
between graphene and Si/Si02 substrates. We mod- 
eled the light modulation by multiple reflections, show- 
ing that: i) the contrast can be tailored by adjusting the 
Si02 thickness. Without oxide, no modulation is pos- 
sible; ii) the light modulation strongly depends on the 
graphite thickness. For few layers (< 6) the samples be- 
have as a superposition of single sheets. For thicker sam- 
ples, both amplitude and phase change with thickness. 
Thus, Rayleigh spectroscopy provides a simple and quick 
way to map graphene layers on a substrate. It can also 
be combined with Raman scattering, which is capable of 
structural identification. 
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